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POPULAR SUMMARY 
HJZ, August 15, 2001 


ICESat’s Laser Measurements of Polar Ice, Atmosphere, Ocean, and Land 


The Ice, Cloud and Land Elevation Satellite (ICESat) mission will measure changes in elevation 
of the Greenland and Antarctic ice sheets as part of NASA’s Earth Observing System (EOS) of 
satellites. Time-series of elevation changes will enable determination of the present-day mass 
balance of the ice sheets, study of associations between observed ice changes and polar climate, 
and estimation of the present and future contributions of the ice sheets to global sea level rise. 
Other scientific objectives of ICESat include: global measurements of cloud heights and the 
vertical structure of clouds and aerosols; precise measurements of land topography and 
vegetation canopy heights; and measurements of sea ice roughness, sea ice thickness, ocean 
surface elevations, and surface reflectivity. The Geoscience Laser Altimeter System (GLAS) on 
ICESat has a 1064 nm laser channel for surface altimetry and dense cloud heights and a 532 nm 
lidar channel for the vertical distribution of clouds and aerosols. The accuracy of surface 
ranging is 10 cm, averaged over 60 m diameter laser footprints spaced at 172 m along-track. 

The orbital altitude will be around 600 km at an inclination of 94° with a 183-day repeat pattern. 
The on-board GPS receiver will enable radial orbit determinations to better than 5 cm, and star- 
trackers will enable footprints to be located to 6 m horizontally. The spacecraft attitude will be 
controlled to point the laser beam to within ± 35 m of reference surface tracks at high latitudes. 
ICESat is designed to operate for 3 to 5 years and should be followed by successive missions to 
measure ice changes for at least 1 5 years. 
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Abstract 

The Ice, Cloud and Land Elevation Satellite (ICESat) mission will measure changes in elevation 
of the Greenland and Antarctic ice sheets as part of NASA’s Earth Observing System (EOS) of 
satellites. Time-series of elevation changes will enable determination of the present-day mass 
balance of the ice sheets, study of associations between observed ice changes and polar climate, 
and estimation of the present and future contributions of the ice sheets to global sea level rise. 
Other scientific objectives of ICESat include: global measurements of cloud heights and the 
vertical structure of clouds and aerosols; precise measurements of land topography and 
vegetation canopy heights; and measurements of sea ice roughness, sea ice thickness, ocean 
surface elevations, and surface reflectivity. The Geoscience Laser Altimeter System (GLAS) on 
ICESat has a 1064 nm laser channel for surface altimetry and dense cloud heights and a 532 nm 
lidar channel for the vertical distribution of clouds and aerosols. The accuracy of surface ranging 
is 10 cm, averaged over 60 m diameter laser footprints spaced at 172 m along-track. The orbital 
altitude will be around 600 km at an inclination of 94° with a 183-day repeat pattern. The on- 
board GPS receiver will enable radial orbit determinations to better than 5 cm, and star-trackers 
will enable footprints to be located to 6 m horizontally. The spacecraft attitude will be controlled 
to point the laser beam to within ± 35 m of reference surface tracks at high latitudes. ICESat is 
designed to operate for 3 to 5 years and should be followed by successive missions to measure 
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1. Introduction 


As noted in the 1992 IPCC Supplement on Scientific Assessment of Climate Change, the largest 
uncertainty about sea level is "rooted in our inadequate understanding of polar ice sheets whose 
response to climate change also affects predictions of sea level rise" [Houghton, et al., 1992], 

The National Academy of Sciences [1990] stated that "possible changes in the mass balance of 
the Greenland and Antarctic ice sheets are fundamental gaps in our understanding and are crucial 
to the quantification and refinement of sea-level forecasts." This uncertainty was further iterated 
by Fitzharris et al., [1996], who noted that "monitoring of key components of the cryosphere 
must continue. The mass balance of the ice sheets of the world is poorly known." Although sea 
level changes in response to input from ice caps and glaciers, as well as thermal expansion of the 
ocean itself, the ice sheets contain the potential for the most significant contribution to sea level 
rise simply because they contain so much ice. The impact of even small rises in sea level has 
been widely studied. Aside from the threat of inundation faced by low-lying coastal areas, 
increased beach erosion that would occur before inundation, for example, is a serious economic 
concern [Yohe and Schlesinger, 1998; Neumann et al., 2000]. 

The primary purpose of ICESat (Figure 1 . 1 and http://icesat.gsfc.nasa.gov/) is the determination 
of inter-annual and long-term changes in polar ice-sheet mass, the causes of changes in mass 
balance (polar precipitation, ice melting, or ice flow acceleration/deceleration), and the impact of 
these changes on global sea level. Changes in ice mass (Figure 1.2) are caused by an imbalance 
between the ice mass inputs (snowfall, condensation, and occasional rainfall) and outputs 
(evaporation, melt runoff, iceberg discharge, and snow drift removal). Conventional methods of 
studying ice sheet mass balance examine the difference between the mass input and output terms, 
but significant errors in these quantities have limited determinations to about ±25% [Warrick et 
al, 1996], which is equivalent to ± 2 mm/yr of sea level change. Measurements of sea level from 
the TOPEX/POSEIDON radar altimeter show a current rise of +2.1 ± 1.3 mm/yr [Nerem et al., 
1997], which is comparable to the rise over the last century measured by tide gauges [Warrick et 
al., 1996], 

The concept of using satellite-altimeters to measure ice elevation changes and determine ice- 
sheet mass balance [Zwally, 1975] is based on the simple relationship of ice surface elevation 
changes to changes in ice thickness and therefore ice mass. Surface elevation changes are 
equivalent to ice thickness changes minus the vertical motion of the bedrock, which is generally 
smaller and can be separately estimated. In addition, short-term changes caused by variations in 
rates of near-surface firn compaction must be accounted for [Arthern and Wingham, 1998 and 
Zwally and Jun, submitted]. Despite the limitations of ocean radar altimeters over the sloping 
surfaces on ice sheet, elevation changes from radar altimeter measurements that are indicative of 
changes in ice mass balance have been reported [e.g. Zwally, 1989; Davis et al., 1 998; and 
Wingham etal., 1998, Zwally and Brenner, 2001]. In 1979, satellite laser altimetry was proposed 
[Campbell et al, 1979] for measurement of ice sheet elevation changes, because of the 


2 


advantages of the smaller laser footprint and the restriction of laser reflections to a small surface 
layer. 

Prediction of sea level rise during the next century requires not only an assessment of the current 
rates of mass gain or loss of the Greenland and Antarctic ice sheets, but information on how the 
ice balance changes with changes in polar precipitation, temperature, cloudiness, and perhaps 
other factors. All of these factors, including the ice sheet surface mass balance processes, 
undergo interannual and decadal variations as well as longer term changes. Therefore, ICESat 
should be followed by successive missions to measure ice changes for at least 15 years. 

In atmospheric science, incomplete knowledge of clouds and aerosols causes one of the main 
uncertainties in modeling and prediction of global climate warming. Climate is especially 
influenced by the vertical structure and horizontal coverage of clouds, which alter the radiative 
fluxes at the top and bottom of the atmosphere and determine the vertical distribution of 
atmospheric heating rates. In particular, direct knowledge of cloud heights has been a limiting 
factor. In addition, while the climatic effect of increasing CO 2 and other greenhouse gases is 
well known, a concurrent increase in anthropogenic aerosols is believed to have an opposite 
cooling effect [Trenberth, 1996]. Aerosols affect the earth’s energy budget and climate by 
scattering and absorbing radiation (direct effect) and altering cloud particle size and number 
density, which ultimately affects the cloud albedo, scattering and absorption properties (indirect 
effect). The exact magnitude of the aerosol cooling is unknown, but is thought to offset a 
substantial fraction of the greenhouse warming. A major reason for the uncertainty has been the 
limited ability to make global observations of tropospheric aerosols, especially over land. 

ICESat will provide continuous global measurements of vertical cloud and aerosol structure and 
optical depth, planetary boundary layer height, and polar tropospheric and stratospheric clouds 
for the first time. The direct and unambiguous lidar data will also be used for validation of the 
cloud and aerosol retrievals from the passive sensors on other EOS satellites. Passive sensors 
view cloud tops very well, but the presence of upper layer clouds limits their ability to 
distinguish multi-level cloud formations and to determine the vertical distribution of clouds. 
Unique aspects of ICESat’ s lidar are the direct measurements of cloud heights, monitoring of 
aerosol distributions over land, and all-year measurements of clouds and aerosols in polar 
regions. In polar regions, improved measurements of clouds is especially important for 
understanding the radiation balance and for improving the atmospheric modeling of 
precipitation, both of which directly affect the surface ice mass balance. 

The land topographic measurements of ICESat will provide a network of high-accuracy profiles 
along the 183-day repeat tracks with a maximum cross-track separation of 15 km at the equator. 
These observations will improve digital elevation models where geodetic control is scarce, 
particularly in the higher-latitude regions. Repeat laser profiling and crossover analysis may be 
used to detect changes in land elevation caused by geologic processes such as soil erosion, 
sediment transport, and magma inflation of volcanoes. Changes in inland water levels can also 
be monitored and vegetation height can be measured in low-relief areas. An important feature of 
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ICESat is the capability to accurately point the spacecraft to off-nadir targets of opportunity, 
such as volcanic aerosol plumes or regions inundated by floods. 

ICESat measurements (Figure 1 .3) will be made by the Geoscience Laser Altimeter System 
(GLAS), which has a 1064 nm laser channel for surface altimetry and dense cloud heights and a 
532 nm lidar channel for the vertical distribution of clouds and aerosols. Laser footprints with 
60 m diameter and 172 m along-track spacing will have a surface-ranging accuracy of 10 cm. 
ICESat’s orbital altitude will be 600 km at an inclination of 94° with a 183-day repeat pattern. 
The on-board GPS system will enable radial orbit determinations to better than 5 cm, and star- 
trackers in GLAS and on the spacecraft will enable footprints to be located to 5 m horizontally. 
Spacecraft pointing of the laser beam can be controlled to within ± 30 m on the surface for 
precise control of repeat tracks over the ice sheets to ± 35 m. 

The ICESat spacecraft and GLAS instrumentation are designed to operate continuously for a 
minimum of 3 years with a goal of 5 years. Full-scale development of GLAS began in May of 
1997 at NASA’s Goddard Space Flight Center, and Ball Aerospace Corporation contracted to 
build the ICESat spacecraft in February 1998. The launch of ICESat on a Boeing Delta-2 rocket 
from Vandenberg, California is planned for December 2001, along with a second payload, 
CATSAT, from the University of New Hampshire. 

2. Ice Sheet Science 

The huge ice sheets of Greenland and Antarctica hold enough fresh water to raise global sea 
level by 80 meters if they melted completely. Although major melting is not imminent, the ice 
sheets change continually in response to natural processes and are expected to respond 
greenhouse-induced climate warming as well. The annual mass exchange between the ice sheets 
and the ocean is about 8 mm/year of global sea level equivalent, so even small changes in this 
rate of mass exchange are significant. Moreover, some parts of the ice sheets may be more 
susceptible to rapid changes. The West Antarctic Ice Sheet alone contains the equivalent of 6-7 
meters of sea level rise and rests on a bed far below sea level, a configuration that may make it 
susceptible to dynamic instabilities [Bindschadler, 1998; Oppenheimer, 1998; Bentley, 1997]. 
The Greenland ice sheet is only about 1/8 as large as the Antarctic ice sheet but may be more 
vulnerable to climate warming because summer melting already occurs over much of its surface. 

Antarctica straddles the South pole and has a dominant influence on its own climate, and on the 
surrounding ocean. Surface melting during summer is mostly restricted to peripheral ice shelves 
and some near coast parts of the ice sheet [Zwally and Fiegles, 1994]. Away from the coast, 
much of Antarctica is a cold desert with very low precipitation rates. Ice drainage is primarily 
through outlet glaciers and ice streams, some of which penetrate deep into the heart of the ice 
sheet moving at speeds of a few hundred meters per year to over 1 km/yr in some cases. Many 
glaciers and ice streams converge into floating ice shelves, which are also fed by snow 
accumulation on their surfaces. Thus, ice loss from the Antarctic ice sheet is primarily through 
the ice shelves by basal melting and iceberg calving at the ice fronts. The ice shelves thin 
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towards their seaward ice fronts, partly by ice creep and partly by basal melting at rates generally 
of a few tens of cm/yr. Although mass changes in the floating ice shelves have no direct effect 
on sea level, but changes in the ice shelves could alter their buttressing effects on the ice streams 
and their rates of discharge from the ice sheets. 

By contrast, the Greenland climate is strongly affected by its proximity to other land masses and 
to the North Atlantic, with the warm Gulf Stream to the south and regions of North Atlantic deep 
water production to the east and west. Summer melting occurs over about 20 % of the ice-sheet 
surface [Abdalati and Steffen, 1997], depending on summer temperatures, with most of the 
resulting melt water flowing into the sea. The average accumulation rate in the accumulation 
zone is about 26 cm water/yr, compared to 1 6 cm water/yr for Antarctica. The zone of net 
ablation, which is about 9 % of the total Greenland ice sheet, has an average annual net mass loss 





water with the ocean is about 25% as large as for Antarctica. Ice-core data from the summit of 
the Greenland ice sheet indicate that Greenland temperatures and accumulation rates can change 
significantly over periods of a few years to decades [Alley et al., 1993]. 

Ice sheets thicken by snow accumulation and thin by snow densification, ice creep and motion, 
and evaporation and melting. Rates of snowfall, evaporation, densification, and melting change 
significantly throughout the year and from one year to the next; ice creep and ice motion exhibit 
only small changes on sub-decadal time scales. Consequently, the ice surface rises and falls in 
response to the seasonal and interannual changes. The magnitude of this short-term elevation 
variability is poorly known, but is estimated to be on the order of tens of cm to a meter over the 
accumulation zones of Greenland and Antarctica. Near the coasts, where snow accumulation and 
melt rates are higher, this short-term variability may increase to a few meters. These short-term 
changes must be monitored for two major reasons: to infer the interannual variability of snow- 


Based on very recent results, the Greenland ice sheet has a negative net balance, contributing 
almost 10% to the total observed sea-level rise [Krabill et al., 2000]. Taken as a whole, the 
higher-elevation parts of the ice sheet are close to balance [Thomas et al., 2000a; Davis et al, 
2000], but there is significant near-coastal thinning [Krabill et al., 2000; Abdalati et al., in press], 
at rates up to 10 m/yr on one major outlet glacier [Thomas et al., 2000b]. The results at higher 
elevations are from three independent analyses, using satellite radar-altimeter data (only south of 
72° N), aircraft laser-altimeter data, and comparison of total ice discharge with total surface 
accumulation. However, the assessment of near-coastal thinning was possible only by the 
airborne laser-altimeter measurements. While the airborne laser system provides accurate 
change measurements in the lower-elevation margins of the ice sheets (where radar altimetry is 
limited by the surface slopes and undulations), the airborne measurements are limited in their 
spatial and temporal coverage. Usually survey flights are made only in the late-spring/early 
summer at intervals of about 5 years. Therefore, the measured changes are necessarily 
characteristic of summer-to-summer trends and may be biased by undetected interannual 
variability. For example, Bromwich et al. [in press] shows increases in precipitation over the 
interior of the ice sheet north of 70° N in the mid to late 1990’s that might not have been 
captured in the airborne-measured elevation changes. This highlights the importance of 
monitoring future behavior of Greenland regions with ICESat data, and of making similar 
surveys over the Antarctic ice sheet where aircraft surveys are not practical. 

The objective of ICESat is to measure profiles of ice sheet surface elevation with sufficient 
accuracy, spatial density, and temporal coverage so that inter-annual and long-term elevation 
changes can be derived with an accuracy of <1.5 cm/yr for spatial-averages of measurements 
over areas of 100 km x 100 km on the ice sheets. The spatial scale is required to resolve details 
of the basin scale mass redistribution and the effects of changes in major ice streams and outlet 
glaciers. The vertical resolution is required to resolve changes that are about 10% of the 
accumulation rate as a typical value. In regions where the accumulation is higher than average, a 
lesser resolution would be acceptable, whereas in some interior low accumulation regions a 
larger spatial averaging will be used to achieve a better than 1.5 cm/yr accuracy. 

The elevation profiles will consist of the mean elevations over the 60 m laser footprints, with an 
elevation accuracy of 15 cm, spaced at 172 meters along-track as described in other sections. 

The pattern of ICESat ground tracks over the margin of west-central Greenland is shown in 
Figure 2.1. Surface elevation changes are usually calculated at crossovers between successive 
ascending and descending passes over the same location (a crossover point), by interpolation 
between successive footprints. Analysis of the error caused by interpolation has been estimated 
with aircraft laser data to be less than several cm over most of the ice sheets where the surface 
undulations between footprints are small. Spatial averages of elevation changes are obtained by 
averaging the elevation differences at a set of crossovers in a chosen area, typically 50 km by 50 
km or larger, taking into account the varying time intervals between the measurements. For 
random errors, approximately 100 crossovers differences are required to reduce the error on the 
mean elevation change to 1 .5 cm. Over most of the ice sheets, the number of available 
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crossovers should be sufficient to achieve < 1.5 cm/yr over 100 km by 100 km areas, using the 
conservative assumption of 50% loss due to cloud cover (i.e. 75% crossover loss). Furthermore, 
in addition to differences at crossovers, considerably more elevation differences will be obtained 
from analysis of the precise-repeat track profiles that will be maintained to within ± 35 m cross- 
track. 

Time-series of surface elevations, H(t), having sufficient resolution to show seasonal changes are 
created by the sequence of average crossover differences between the first 90-day interval and 
each of the successive 90-day intervals, combined with the sequence from the second interval 
crossed with each successive interval, and so forth for the sequences for the third and greater 
intervals [Zwally and Brenner, 2001]. As the length of the measurement period increases, the 
number of crossovers for analysis increases in proportion to time squared. For an average of M 
crossovers in each pair for N time intervals, the total number of crossovers in the series for the 
first sequence from crossovers with the first interval is N x M. However, the total number of 
crossovers in the combined H(t) series for all sequence is M x £[ N (n) = M x N(N+l)/2 , which 
is (N+l)/2 times larger. In addition to the traditional crossover differences, elevation differences 
from the precise-repeat track profiles can be included in the time series in a similar manner. The 
H(t) time series is then analyzed to infer the character of the changes, for example, a seasonal 
cycle imposed on a linear trend. 

H(t) time series constructed from crossovers from 7 years of ERS- 1 and ERS-2 data at two 
locations in West Antarctica, are shown in Figure 2.2. In each case, crossovers within a 100 km 
radius and ± 250 m elevation of the central point are included. The data are corrected for an 
unexpected inter-satellite bias that was determined by Brenner et al. [2000] from analysis of 
crossover differences acquired during the 12 months of overlapping operation of ERS-1 and 
ERS-2. The bias was found to vary with received backscatter power, which depends on the 
surface slope and therefore varies with elevation. In this region, the bias lowers the ERS-2 
elevations relative to ERS-1 by 17.4 cm at 500 m elevation and 1 1 cm at 2000 m. 

The first H(t) series, which is centered at Byrd Station (80° S, 120° W) at 1526 m, shows a linear 
increase of 1 1.9 cm/yr and a peak-to-peak seasonal amplitude of 67 cm. The second location 
(76° S, 105° W) at 892 m elevation about 100 km inland of the Thwaites Glacier tongue has a 
significant thinning of -20.5 cm/yr, with only a small and irregular seasonal amplitude of 13 cm. 
For the ice sheet summit (72.57° N, 38.45° W) in Greenland, a fitted seasonal amplitude of 20 
cm with more interannual variability than at Byrd has been modeled [Zwally and Jun, submitted] 
using a temperature-dependent rate of fim compaction and the observed seasonal and interannual 
variations in temperature from automatic-weather-station data [K. Steffen, Pers. Com.]. Most of 
the annual compaction of the upper layers of fim occurs during the warmer summer months, as 
illustrated by the March minimum in the fitted seasonal amplitudes in Figure 2.2. 

Previously, Wingham et al. [1998] reported a -1 1.7 ± 1.0 cm/yr thinning trend averaged over the 
Pine Island Glacier and Thwaites Glacier drainage basins for the period 1992-1996, also from 
ERS data. This region is of particular interest, because Pine Island and Thwaites Glaciers drain a 
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large part of the West Antarctic ice sheet and may be undergoing significant changes [Vaughan 
et al, 2001]. For example, recent analysis of SAR interferometry showed the position of the 
hinge line of the Pine Island glacier retreated 1.2 km/yr in the same period as the altimeter 
measurements of elevation change [Rignot, 1998]. In the region from 1 10° to 130° W, the 
elevation changes range from about - 30 cm/yr to + 25 cm/yr, as shown in the map of dH/dt in 
Figure 2.3. The dH/dt increases are generally at the higher elevations and along the ridge 
between the Pine Island and the Ross-Ice Shelf drainage basins. 

During 3 to 5 years of ICESat-1 operation, an estimate of the overall ice sheet mass balance and 
sea level contribution during the period of measurements will be obtained. One estimate of the 
accuracy of a linear trend deduced from ICESat data over 5 years for all of Antarctica is ~ 7 
mm/yr of water equivalent, including errors in the estimates of fim compaction and post glacial 
rebound [Wahr et al., 2000]. Their estimate of the accuracy for deducing long-term trends was 9 
mm/yr, obtained by including the difference between a typical 5-year and century-scale trends 
from simulations of variations in Antarctic precipitation. 

Currently, correction for post-glacial rebound is made using coupled models of ice sheet 
evolution and the visco-elastic response of the solid Earth [Huybrechts and Le Meur, 1999]. In 
the future, however, estimates of the basal motion can be improved [Wahr et al., 2000] through 
the combined analysis of data from ICESat and the Gravity Recovery and Climate Experiment 
(GRACE) missions, which are expected to fly concurrently. Gravity measurements are sensitive 
to changes in mass under the satellite, including mass due to crustal uplift and mantle inflow and 
the mass of snow and ice in the ice sheet. For a given mass change, the surface height changes 
would be in the ratio of 1 to 3.5 to 10, respectively, for isostatic uplift, changes in solid ice, and 
changes in snow accumulation at the surface because of their different mass densities. 

Therefore, measurements of both change in height and change in gravity make it possible to 
improve the post-glacial rebound estimates, with an accuracy that improves with the length of 
concurrent measurements. 

The importance of continued ice-elevation change measurements after the first ICESat is 
illustrated by the variability in the area of Greenland surface melt observed from satellite passive 
microwave data over 17-years and its correlation with temperature Figure 2.4. Whereas the 
passive microwave only provides the area of surface melting, ICESat’s elevation changes will 
provide the volume of melt water produced from the marginal ablation zones. While rates of ice 
flow in the ice sheets change only slightly on decadal time scales, the surface mass balance and 
mass exchanges with the atmosphere and ocean change instantaneously with changes in 
precipitation and melting. For the purpose of understanding and predicting changes in ice mass 
and sea level, observation of the relationships between changes in ice melting and accumulation 
and concurrent changes in climate over decadal time scales is essential. ICESat data on seasonal 
and inter-annual changes in surface elevation will provide unique information for validation of 
atmospheric models of precipitation-evaporation in polar regions, as well as energy-balance 
models of surface melting and snow accumulation on the ice sheets. 
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3. Atmospheric Science 


For several decades, the value of a satellite lidar to directly measure the height distribution of 
clouds and aerosol in the atmosphere has been considered [Curran, 1989]. Current satellite and 
other data provide considerable information on cloud and aerosol distributions [e.g. Rossow and 
Cairns, 1995]. However, some critical parameters can only be obtained by active optical 
profiling. Specifically, passive techniques can not resolve the height profile of aerosols, which is 
required for modeling aerosol transport and height-resolved radiative heating/cooling effects. 

For transmissive and multi-layer clouds, passive techniques are limited in accuracy. The ICESat 
GLAS active lidar profiling will supplement passive cloud sensing with data to vertically resolve 
the heating and cooling effects of global cloud cover and aerosols. 

The concept of the GLAS atmospheric lidar measurements as a combined measurement with 
laser surface altimeter is based on the fundamental compatibility of basic instrument 
requirements. Although atmospheric signals are typically much smaller than for the surface 
return signals, the vertical resolution requirement is much less; 75 m is sufficient for almost all 
requirements. As is well known, the signal to noise of measurements increases with decreasing 
bandwidth or resolution. Thus there is compatibility in the laser radar, power-aperture product 
requirements between precision surface altimetry and basic profiling of clouds and aerosol in the 
atmosphere. The cost of a space borne mission to obtain both observations is significantly less 
than for two separate experiments. 

3.1. Atmospheric Measurement Objectives 

The atmospheric science objectives for cloud heights and vertical profiles, polar clouds, aerosol 
profiles, planetary-boundary-layer heights, and cloud filtering for surface altimetry are described 
here. The GLAS measurement of cloud heights and vertical profiles are described in detail in 
Spinhime and Palm [1996]. The along-track laser measurements provide an image of the 
vertical structure, or morphology of cloud systems as shown in Figure 3.1 using aircraft data 



tropopause. Although the optical thickness of this layer is small, the net radiative impact over 
the entire region is estimated to be significant. Reliable detection and appropriate sampling is 
not possible by passive means. The GLAS lidar will provide a very sensitive measure of the 
presence, height and thickness of tenuous cirrus layers. 

A major goal is to observe multiple-layered and broken-cloud cover, as well as the height of the 
cloud tops. For a large fraction of cloud cover, GLAS will detect the cloud-base height. Cloud 
effects on surface warming are largely related to the cloud component of the downward 
atmospheric flux. A principal determining factor of the cloud flux component is the distribution 
of cloud-base heights. For clouds forming in the PBL, the cloud base indicates the lifting 
condensation level (LCL). The LCL is a significant variable for boundary layer dynamics. From 
the LCL, the near surface moisture over the oceans can be accurately derived [Palm, et al., 

1998 ], 

Polar cloud observations will have special value, because passive techniques are limited by 
darkness part of the year and cold temperatures and bright backgrounds (snow and ice) the rest of 
the year. GLAS measurements will unambiguously define cloud type and fraction, which define 
the net radiation balance of the Arctic regions and strongly affect atmospheric dynamics. In the 
cold seasons in the Arctic, hazes of small ice crystals, or diamond dust, occur in the lower 
troposphere greatly reducing visibility. The ice crystals are an important factor in the radiation 
transfer, which determines the vertical temperature and humidity structure of the winter Arctic 
atmosphere. For studies of ice sheet mass balance, cloudless ice-crystal precipitation is believed 
to be a dominant source of precipitation in large areas of Antarctica. GLAS active laser sensing 
will give the first spatial and temporal distribution of clouds and diamond dust for polar 
precipitation. 

Aerosol profiles provided by the GLAS instrument will uniquely define the vertical aerosol 
structure throughout the troposphere over land and ocean. Vertically-elevated layers of 
particulates are transported over long ranges and have been linked to major climate and 
atmospheric-chemistry impacts. Also, transported aerosols are thought to influence air quality 
and ocean nutrients. Ground-based lidar measurements have long been a standard for monitoring 
volcanic aerosol layers in the stratosphere. The ICESat measurements will give global coverage 
of episodic aerosol events such as volcanic emissions, biomass burning, ablated desert soils, 
continental particulates, and arctic haze. When coupled with wind direction and speed, the 
vertical distribution of aerosol trace materials provides information on aerosol mass transport. 
Long-range transport of trace gases and aerosols is a dominating factor in the global chemical 
balance, for which accurate aerosol vertical-profile information is required. 

The planetary boundary layer height (PLB) is a basic parameter linking the surface to 
atmosphere dynamics, especially over the world oceans. However, there is currently no way to 
obtain accurate PBL height measurements on a global basis. GLAS observations will determine 
the heights of cloud-capped boundary layers with high accuracy. In many areas, GLAS 
observations will also determine the height of the mixed layer from the aerosol scattering 
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In order to fulfill the science requirements, GLAS must measure the atmospheric scattered signal 
at the accuracies and spatial scales given in Table 3a. These values are based on a complex 
analysis and extensive experience with aircraft lidar remote sensing [Spinhime and Palm, 1996]. 
The obtainable measurement accuracy and horizontal and vertical resolutions depend upon the 
instrument configuration and the atmospheric scattering cross sections. In essence, both 
horizontal and vertical averaging of data is available to make a given measurement when 
atmospheric cross sections are 10' 4 (m-sr) -1 or smaller. 

Model-derived estimates of the minimum detectable backscatter cross section as a function of the 
horizontal resolution (or number of shots averaged) are shown in Figure 3.2 for nighttime and the 
worst expected daytime background conditions. During nighttime, GLAS is expected to resolve 
thin (sub-visible) cirrus layers and tenuous aerosol layers (5x1 O' 8 (m-sr) ') at a horizontal 
resolution of about 50 km. Normal cirrus and denser aerosol layers (10' 6 (m-sr) _l ) will be 
detectable at about the 2 km resolution. Clouds with backscatter cross sections greater than about 
10' 5 (m-sr) 1 will be detected on a shot to shot basis (172 m horizontal resolution). During 
daytime, these resolutions will vary according to the background intensity, but will generally be 
about 5 to 10 times less than the nighttime resolution. For horizontal averaging intervals less than 
1000 km, the absolute minimum detectable backscatter cross section is about 10' 7 and 10 9 (m- 
sr)' 1 during day and night, respectively. 


4. Land and Vegetation Science 

The Earth’s surface consists of a complex mosaic of land forms and vegetative cover resulting in 
large variations in elevation, slope, roughness, reflectance, and vegetation height. Present day 
landscapes have been generated by the integrated effects of a diverse set of lithospheric, 
cryospheric, hydrospheric, ecologic, atmospheric, and anthroprogenic processes. Documentation 
of these landscape properties, and the associated dynamics, is a first step in understanding the 
interplay between the formative processes. Characterization of the landscape is also necessary 
for physical models that depend on interactions between the land and other components of the 
Earth system, such as predictive models of global climate change. Adequate knowledge of these 
characteristics is also essential for proper utilization of the land as a resource, as well as for 
prediction and mitigation of natural hazards (e.g., earthquakes, volcanic eruptions, landslides, 
floods, wild fires). 

Current physical descriptions of land properties on a global basis are in many respects 
inadequate. The best publicly available global representations of the Earth’s topography are 
Digital Elevation Models (DEMs) with a spatial resolution of 1 km [GLOBE Task Team and 
others, 1999], which are too coarse for the needs of most physical studies of land process. 
Moreover, these DEMs are amalgams of many disparate topographic sources, often with 
inconsistent or poorly known accuracies and a bewildering collection of horizontal and vertical 
datums. Global characterizations of land cover [e.g., Hansen et al., 2000; Loveland et al., 2000] 
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and vegetation properties [e.g., DeFries et al., 2000] have recently been produced at moderate 
spatial resolution using satellite remote sensing data. Comparisons of alternative global land 
cover products [Hansen and Reed, 2000] reveal considerable discrepancies between the products, 
and there is currently no satisfactory method for validating global maps of land cover or 
vegetation properties. 

ICESat will repeatedly sample global land topography and vegetation heights to characterize 
landscape properties and monitor changes in dynamic landscapes. Of particular significance is 
the ability to re-profile ground tracks very accurately by off-nadir pointing to specific targets of 
interest. In mid-latitudes (N to 50°), off-nadir pointing will normally not be used to follow a 
prescribed ground track, so tracks on repeat cycles will vary within about 1 km. However, off- 
nadir pointing up to 5° to profile specific targets might be done several times a day. At the 
equator, a selected location could be targeted up to 12 times during the 183-day orbit cycle, and 
more often at higher latitudes. After launch, potential operational limitations on pointing 
commands will be better defined. Potential applications include long-term monitoring of water 
levels of selected rivers and lakes, retreat or surge of alpine glaciers, soil erosion, changes in 
snow depth, inflation or deflation of volcanoes, and regional response of the lithosphere to ocean 
tides, changes in ice loads, and tectonic stresses. Also, the effects of natural or man-made 
catastrophes on topography and vegetation can be targeted. Alternatively, ICESat’ s precise 
pointing capability can be used to target slightly-offset ground tracks within a selected region to 
achieve a more finely resolved grid of profiles for more detailed local mapping. 

ICESat elevation profiles will provide independent high-accuracy data to calibrate and validate 
topographic and land cover products generated by other means. In particular, high-resolution 
DEMs and land cover products are being derived from optical imaging and synthetic aperture 
radar (SAR) systems. For example, the recent Shuttle Radar Topography Mission (SRTM) will 
yield a 30 m horizontal-resolution DEM between 60° N and 57° S latitudes from interferometric 
SAR, and the ASTER instrument on Terra yields local DEMs of similar resolution by stereo- 
photogrammetry. Typical vertical resolution is a few meters. However, these DEMs have 
inherent systematic and random errors. The ICESat elevation profiles, defined in a consistent. 
Earth-centered reference frame, can be used to evaluate these DEMs and establish global 
geodetic control, as demonstrated using profiles from the Shuttle Laser Altimeter experiment 
[Sun and Ranson, 1997; Garvin et al., 1998; Harding et al., 1999; and Sun et al., 2000]. 

Derivation of vegetation height data, and topographic profiles of the ground surface where it is 
covered by vegetation, depends on interpretation of the GLAS altimetry channel waveform 
record for each laser pulse. ICESat waveform processing (Section 7) will decompose the 
waveforms into vertically distinct components by fitting multiple Gaussian distributions, 
following methods developed using waveforms acquired by the Shuttle Laser Altimeter 
[Carabajal et al., 1999] and the airborne Laser Vegetation Imaging Sensor [Hofton et al., 2000a]. 
Interpretation of the resulting Gaussian fits will depend on land cover properties and local 
surface relief. For footprints lacking vegetation cover or cultural features, the overall centroid of 
the Gaussians fit is taken to be the mean surface elevation within the footprint. In the case of 
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multiple returns, the individual Gaussians are assumed to define the elevation of discrete 
surfaces: for example, upper and lower flat-lying areas separated by a cliff. The width of each 
Gaussian contains information on the relief of the corresponding surface due to the combined 
effects of slope and roughness. For footprints containing vegetation and/or cultural features, 
surface relief effects are combined with those of the height distribution of canopy components 
(living or dead foliage and woody tissue) and cultural features. Where the surface relief is flat or 
very muted and sufficient laser energy penetrates the canopy to the underlying ground, a lowest 
return is acquired in the waveform and is interpreted as originating from the ground surface. 
Higher returns are interpreted to correspond to vegetation layers and/or the tops of cultural 
features. The very first return corresponds to the upper-most detected canopy surface or highest 
cultural feature. Where surface relief is pronounced the ground, vegetation, and cultural features 
become unresolved in the waveform because their height distributions intermix. 

The use of laser waveforms to derive vegetation height information and ground topography 
beneath vegetation has focused on contiguous footprints with diameters in the range 8 to 25 m 
acquired by airborne systems [Blair, et al, 1999; Harding et al., 2000], Hofton et al. [2000b] 
describe topographic mapping of the Long Valley caldera in California and Weishampel et al. 
[2000] report on canopy vertical structure mapping of a tropical rainforest in Costa Rica. Laser 
waveforms have been used to estimate ecologically important forest stand attributes including 
stem density and above-ground biomass [Lefsky et al., 1999a, 1999b; Means et al., 1999], and 
the vertical distribution of plant area [Harding et al.. In press] and light transmission [Parker et 
al., In press]. Unlike the smaller, contiguous footprints of the upcoming Vegetation Canopy 
Lidar mission [Dubayah et al., 1997], the large, non-contiguous GLAS footprints are not 
optimized for retrieval of vegetation height data. Nonetheless, the GLAS waveforms will 
provide important information on vegetation height, especially in areas of low relief where the 
ground and canopy returns can be decomposed. Harding et al. [in press] computed average 
canopy height profiles for areas nominally 50 x 50 m in size, nearly the size of a GLAS footprint. 
These large-area height profiles reproducibly revealed variations in ecologically important 
canopy attributes, including maximum height, and the height, depth and relative plant area of 
overstory, midstory, and understory layers. In sloped or topographically rough areas, ICESat 
will provide information on the total within-footprint relief due to the combined vertical 
distributions of the ground, vegetation, and cultural features. 


5. Sea Ice and Ocean Science 

The polar oceans are covered by a seasonally varying layer of sea ice, varying in thickness from 
centimeters to a few meters [e.g. Gloersen et al., 1992]. Sea ice affects both the overlying 
atmosphere and the underlying oceans in several ways. It is a strong insulator, limiting heat 
exchange between ocean and atmosphere, it modulates the exchange of momentum between 
atmosphere and ocean, and with its high albedo it strongly affects the absorption of radiant 
energy by the Earth. During formation, sea ice rejects salt, and it produces fresh water when 
melting, thus affecting the salinity structure of the ocean with important ramifications for deep 
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convection, bottom-water formation, and blooms of ocean biological productivity associated 
with the ice edge in spring. In addition to these characteristics that influence weather and climate 
in ways that are still poorly understood, sea ice obstructs shipping, modifies submarine acoustics, 
and provides an environment essential to the survival of a wide variety of polar animals. 

Sea ice is a complex material consisting of an ice matrix with inclusions of air, brine, solid salt 
and contaminants. It varies, spatially and temporally in thickness, composition, snow cover, 
wetness, and surface roughness. During the winter, the surface is generally covered by snow, 
which in summer may melt to form slush and melt ponds. Particularly in the Antarctic, where 
ice is thinner, heavy snow loads can depress the ice sufficiently to permit flooding of seawater 
above the ice/snow interface. Since 1972, satellite remote sensing systems have been the major 
tools for mapping and monitoring the sea ice, with emphasis on microwave sensors, both active 
and passive [e.g. Carsey, 1992; Comiso, 1995; Perovich, 1996]. However, sea ice thickness and 
roughness are important parameters for which observation from space has been elusive. 

Most sea ice can be regarded as a nearly-horizontal rough surface with ridges and hummocks, 
which are formed during convergence of the ice pack, and leads with open water or thin ice, 
which are formed during divergence of the ice pack. The roughness is an indication of the 
history of the ice, and strongly influences its drag coefficient and hence its response to winds in 
terms of turbulent energy exchanges and drag-induced flow dynamics. The sea-ice elevations 
will help improve knowledge of the geoid in polar regions. The reflectance is indicative of sea- 
ice albedo, which also affects energy exchange. The surface elevation of flat regions is slightly 
higher than local sea-surface elevation, depending on the ice thickness and snow cover. Sea-ice 
covered ocean also contains icebergs, with surface elevations of tens to several hundred meters. 
The major sea-ice parameters to be determined from GLAS ranges and waveforms shapes 
include surface elevation, surface roughness, and reflectivity. In addition, the freeboard of 
icebergs will be measured wherever ICESat orbits pass over icebergs. Recent work with satellite 
radar-altimeter data [Peacock et al, 1998] indicates the possibility of estimating sea-ice 
freeboard, a proxy indicator of sea ice thickness by comparing elevations of open-water leads 
with those of the intervening ice. Data from GLAS could be better suited to this application, 
because of their small footprint and unambiguous measurement of the snow-ice surface over all 
types of sea ice. Even if identification of open-water waveforms is difficult, it might be possible 
to identify leads by the abrupt change in elevation between lead and sea ice. If successful, this 
capability would represent a major enhancement to our ability to monitor sea ice thickness and 
climatic changes. 

ICESat will spend much of its time over the ocean acquiring information on sea-surface 
characteristics. The shape of the GLAS return-pulse waveform will be determined primarily by 
the surface-height distribution within the footprint, which is small enough to be affected by 
individual large waves. Over distances of cm to a few hundred meters, the sea surface is 
roughened by waves and ocean swell, but over distances of many km the sea surface is almost 
flat. Nevertheless, surface slopes and long-wavelength undulations are present, caused by 
variations in Earth's gravity field, associated for instance with sea mounts, ocean currents, and 
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variations in atmospheric pressure and seawater density. Satellite radar altimeters have shown 
remarkable success in measuring sea-surface elevation and significant waveheights [e.g. L. Fu et 
al, 1994]. Because of its large footprint (> 1 .5 km), a satellite radar altimeter averages the effect 
of the small-scale roughness in the information in the composite return pulse along 300 to 600 m, 
from which estimates of surface elevation and roughness are derived. Consequently, the resulting 
sea-surface elevations can be used to study the longer-wavelength variability, and estimates of 
surface roughness are a statistical indication of the wave height. 

Each GLAS return pulse from the ocean surface will provide estimates of average surface 
elevation and of surface-height distribution within the corresponding 60-meter footprint, at a rate 
of 40/second. However, these footprints will, in general, cover less than one of the longer ocean 
waves, and it will be necessary to include information from many consecutive footprints to infer 
sea-surface elevation and wave height. Nevertheless, we expect GLAS to complement satellite 
radar-altimeter data, particularly near coastlines and islands, where the large radar beam-limited 
footprints include the effects of coastal topography. Moreover, details of the surface-height 
distribution within the laser footprints may help correct for electro-magnetic bias errors in the 
radar data caused by asymmetries in the wave-height distribution. 


6. GLAS Instrument 

The GLAS instrument will measure the distance to the Earth’s surface with 1064 nm laser pulses 
and the vertical distributions of clouds and aerosols in the atmospheric with both the 1064-nm 
pulses and 532 nm pulses. The measurements include the strengths of the received signals and 
the height distributions of the reflecting surfaces. GLAS [Abshire et al., 2000] includes three 
identical diode-pumped Q-switched ND:YAG laser transmitters [Afzal et al., 2000], a 100 cm 
diameter Beryllium receiver telescope, solid state detectors for 1064 and 532 nm signals, a 
subsystem to measure the pointing angle of each laser pulse, and waveform digitizers to record 
the laser backscatter signals at both 1064 and 532 nm. The GLAS components (Figure 6.1) are 
mounted on an L-shaped optical bench and the thermal radiators are mounted at right angles on 
the ends of the benches. The instrument specifications are summarized in Table 6. 1 . 

GLAS has 3 identical Nd:YAG laser transmitters to meet a 3-5 year lifetime requirement. One 
laser will operate continuously, and alternate lasers will be turned on and optically selected via 
flip mirror assemblies when needed. The lasers are passively Q-switched, diode-pumped, and 
conductively-cooled. Approximately 6-nsec wide pulses are emitted in a TEMoo beam. The 
laser system uses a - 1 mJ energy Q-switched master oscillator to establish the pulse length and 
spatial beam quality. Two double-pass, zig-zag slab laser amplifiers increase the pulse energy. 
The total energy output from amplifier at 1064 nm is 75 mJ, of which about 50 mJ is transmitted 
and 25 mJ is directed to a non-linear crystal to produce 35 mJ of transmitted energy at 532 nm. 
The theoretical estimate of the GLAS ranging error is well under 10 cm for clear-sky 
atmospheric transmissions of 50% and surface slopes under 5° (Figure 6.2). The transmitted 
pulse has a 0. 1 10 urad divergence, which will illuminate a spot on the Earth’s surface with a 
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diameter of about 60 m. The laser pulses at 40 Hz repetition rate, producing a separation of 
illuminated footprints on the Earth’s surface of 172 m. Ten-second time markers from the GPS 
receiver provide the long-term standard for the GLAS clock frequency. 

The laser backscatter signals are collected by the Cassegrain-design telescope with a 700 urad 
field of view. The 1064-nm receiver has an 800 pm optical bandpass and two (redundant) Si 
APD detectors with -160 MHz bandwidths. The second detector can be selected if needed via a 
flip mirror. The 1064 nm receiver digitizes both the transmitted and backscatter signals at 1 GHz 
for analysis of the surface returns (Figure 6.3). The 1064 nm signal is also digitized at 2 MHz 
for detection of dense clouds and aerosols at 75-m vertical and 172-m horizontal resolutions. 

The backscatter measurements at 532 nm are more sensitive to weak scattering layers in the 
atmospheric, due to the utilization of photon-counting techniques and along-track signal 
integration. Eight sensitive photon detectors and a narrow 30-pm etalon filter enable the 
detection of optically-thin aerosols during both day and night conditions. Accumulation of 20 
consecutive 532-nm backscatter profiles at 40 Hz will produce an average vertical profile at a 2 
Hz rate with an along-track resolution of 3.44 km. 

The pointing angle of each laser beam is measured relative to inertial space to better than 10 
urad accuracy by the GLAS stellar reference system (SRS), as described by Sirota et al. [2000]. 
The SRS measures the inertial orientation of the rigid GLAS optical bench using a zenith- 
viewing star camera measurement at 10 Hz and a gyroscope for noise filtering. A laser reference 
system (LRS) measures the orientation of each 40-Hz laser beam with respect to the optical 
bench. The LRS extracts a small fraction of the outgoing laser beam into the laser reference 
camera (LRC) with two highly-stable cube corners. A small laser-profiling array digitizes the 
far-field pattern of the laser beam at 40 Hz to determine its orientation to the LRS. A CCD 
camera on a telescope of the LRS digitizes zenith-sky images to determine the inertial orientation 
of the LRS at the same 10-Hz rate as the star camera, but with a higher resolution and a narrower 
field-of-view. Common stars will appear in the LRC and the star camera at an average interval 
of 5-10 minutes for cross-referencing. In addition to the star tracker on the GLAS optical bench, 
the bus star trackers are available as a backup. Although, the mounting of the GLAS instrument 
to the spacecraft will distort in response to thermal variations orbit, the distortions will be 
monitored in-flight and thermally modeled. 

For each laser pulse, about 4,500,000 1-ns digitized waveform samples of the 1064 signal with 
range resolution of 15 cm per sample will be temporarily stored on-board. The samples will be 
searched for the transmitted pulse and for the surface-reflected signals in order to select the 
relevant samples to send to the ground. Start and end points of a search window for the surface 
return (up to 1 1 km in range) will be set for each pulse using a global DEM to predict the surface 
range. The surface finding algorithm passes the digitized waveform samples through six low- 
pass digital filters to maximize the probability of finding echoes from sloped or rough terrain and 
to minimize the probability of selecting cloud echoes. A separate threshold is set for each of the 
filtered waveforms to distinguish the signal level from noise, based on a noise estimate from the 
samples a 1 km region beyond the end of the range window (i.e. below the predicted surface). 
Because the surface return is expected to be the last maximum in the range window, the 
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algorithm searches for the surface echo backward in time from the end of the search window 
towards the start. The echo pulse, defined by the first threshold crossing (from below the 
threshold to above) to the next threshold crossing (from above the threshold to below) is selected 
for each filter for which such a pulse exists. Once the filter has been selected, it is used to 
determine the start of 1000 selected samples. Normally the first 544 or these samples over ice 
sheets and land (8 1.6 m range) or the first 200 samples over sea ice and oceans (30 m range) will 
be further selected for transmission at the full 15-cm resolution. If desired to enlarge the 
measured range span to cover tall trees or areas of large topographic relief, the samples can be 
averaged by factors of 2, 4, or 8. 


Table 6. 1 . GLAS Instrument S 

pecifications 

Laser Type 

ND:YAG slab, 3 stage 
Q-switched, Diode-pumped 

Number of lasers 

3 (one on at a time) 

Laser firing rate 

40 pps 

Laser pulse width 

6 nsec 

Laser Energy (nominal) 

75 mJ, 1064 nm (total) 

50 mJ, 1064 nm (transmitted) 
35 mJ, 532 nm 

Laser Divergence angle 
( 1/e 2 pts) 

1 1 0 urad 

Telescope diameter 

100 cm 

1064 nm detector 

Si APD - analog (1 + 1 backup) 

1064 surface digitizer 
resolution 

1 nsec 

1064 nm cloud digitizer 
resolution 

500 nsec 

532 nm detector 

Si APD - Geiger 
(8 each) 

532 nm aerosol digitizer 
resolution 

500 nsec 

Beam inertial orientation 

< 10 uradian 

Mass 

330 kg 

Power 

3 1 0 W average 

Instr. Duty cycle 

100% 

Data rate 

~ 500 kbps (uncompressed) 

Physical size 

~ 1 10 x 140 x 1 10 cm 

Thermal control 

Radiators with variable 
conductance heat pipes 


The atmospheric backscatter profiles at 1064 nm are also recorded for every laser pulse to 30 
km above the surface. The backscatter profile is electronically filter to 2 MHz (75 m) resolution. 
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Two approaches are used to limit the noise from daytime background light. The 532 nm receiver 
has a narrower 170 urad field of view, and it incorporates a movable mirror assembly to keep the 
FOV centered on the laser signal. A 370-pm wide optical pre-filter and Fabry-Perot etalon also 
reject background light. The etalon has a bandwidth of 30 pm, and its center wavelength is 
thermally tuned via a feedback loop to track the average 532 nm laser wavelength. The 532 nm 
laser backscatter signal that passes through the etalon filter is distributed via beam splitters to 
eight photon counting detectors, which operate in the Geiger mode with ~ 60% photon-counting 
efficiency. 


7. Surface Return Waveform Processing 

The range to the mean surface within the laser footprint is determined from one-half the 
measured time from the transmitted pulse to the time of receipt of the pulse reflected from the 
surface taking into account the spreading of the respective pulses (Figure 7.1) . The transmitted 
pulse has a distribution in time that is approximately Gaussian with a or of about 3 nsec. The 
shape of the return waveform is affected primarily by the transmitted pulse shape, the surface 
height distribution within the footprint, and atmosphere scattering. The surface is characterized 
in terms of roughness and a mean slope across the footprint. (Figure 7.2). Since the digitized 
transmitted pulse is Gaussian, the return shape will also be Gaussian if the surface height 
distribution is Gaussian Gaussian and atmosphere scattering effects are small. Although the 
surface height distribution is not always Guassian, analysis of numerous waveforms acquired by 
the shuttle laser altimeter has shown that most waveforms are approximated very well by a 
Gaussian or as a sum of Gaussians plus a noise bias. 

Over the oceans, sea ice, and most of the ice sheets, the return is expected to be a single 
Gaussian. Over land and more complicated ice sheet regions, there may be multiple distinct 
surfaces within the footprint that will show up as multiple peaks in the return waveform. 
Therefore the modeled waveform is defined as 


w(t) —e exp 

m = 1 


0 — ) 

2C Jn 1 


where N p is the number of peaks found in the waveform, e is the bias (noise level), A m is the 
amplitude, t m is the time, O m is the standard deviation of the m th peak. 

Nonlinear least squares fitting will be used to compute the model parameters in Eq 15. The 
maximum number of peaks (N p ) normally calculated will be six. The range to the surface is the 
time between the centroid of the transmitted pulse and the center of the Gaussian fitted to the last 
peak of the received pulse. Forward scattering from cloud layers in the atmosphere cause some 

19 



photons to take a longer path, and the consequent time delay causes an asymmetric waveform 
with a long tail to the right for the strongest scattering cases. However, since the range is 
calculated to the center of the fitted Gaussian rather then to the centroid of all the data 
surrounding the peak, the effect of forward atmospheric scattering on the calculated range is 
significantly reduced. The surface height distribution representing the combined surface 
roughness and slope is calculated from the standard deviation, Om, of the Guassian for the 
surface peak. Since it is not possible in general to separate slope and roughness, our algorithm 
will calculate both roughness assuming zero slope and slope assuming zero roughness. 


8. Spacecraft and Orbit Characteristics 

The ICESat spacecraft is a version of the commercial RS-2000 spacecraft by Ball Aerospace in 
Boulder, Colorado. The RS-2000 bus, which was designed for acquiring high-resolution (1-m) 
commercial imagery, was first used for Quikscat launched in November 1998. The ICESat bus 
mass is about 670 kg and the instrument mass is about 330 kg, for a total mass of about 1000 kg. 
Orbital average power capability (5 yr end of life) is 630W to provide 330 W for GLAS, 216W 
for the bus, and a 84W margin. The total area of the solar-panel arrays of triple junction cells is 
about 8 m 2 . Although the solar panels rotate, periodic yaw maneuvers of the spacecraft are 
required to reorient panels with respect to the sun. In “sailboat” orientation, the velocity vector 
is parallel to the solar panel axes and in the “airplane” orientation the velocity vector is 
perpendicular. The maximum area-to-mass ratio of the observatory is about 0.01 m 2 /kg. 

The spacecraft carries instrumentation for the near-real time on-board orbit propagator and 
precise post-processing orbit determination. Two redundant Blackjack dual-frequency Global 
Positioning System receivers are connected to separate antennas mounted on the zenith deck. 

The receivers are similar to those on the CHAMP, GRACE and JASON missions. A unique 
feature is the ability to track the GPS signals on both L-band frequencies for precise positioning. 
The hardware and software is capable of tracking 9 or more GPS satellites. Additional orbit 
positioning is provided by a passive laser reflector array mounted on the nadir deck for ranging 
from ground-based laser systems, primarily for validation of the GPS positioning. 

The spacecraft carries two Ball CT-602 star trackers mounted on the zenith deck in orthogonal 
directions. In flight, the star-image data at 10 Hz and the gyro data from the GLAS instrument 
are processed in the spacecraft computer to determine the ICESat attitude autonomously for 
spacecraft attitude control. The agile attitude-control-system on the spacecraft allows near-nadir 
pointing control in normal operation and off-nadir pointing by rolling and/or pitching the 
spacecraft. The accuracy of the spacecraft pointing control near nadir is ± 10 arcsec (lO), which 
is equivalent to ± 30 m on the Earth. The sub-orbital track spacing will be controlled to ±1 km at 
the equator by periodic spacecraft thrusting [Demarest and Schutz, 2000]. Based on predictions 
of solar flux at launch, the frequency of the maneuvers is expected to be about 5 days. At 
latitudes above 50° N and S, the spacecraft will be rolled slightly to point to the reference ground 

tracks defined by the orbit repeat cycles to an accuracy of ± 35 m (10) on the surface, including 
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the errors in pointing the predicted orbit position between pointing-command uploads to the 
spacecraft. Globally, for targets of opportunity, the spacecraft can be commanded to off-nadir 
angles, normally limited to ± 5° or about ± 50 km on the surface. 

The 94° orbit inclination is a compromise between the desire for true polar orbit to observe the 
total Antarctic ice sheet and the need for orbital crossovers for elevation change analysis. A 
crossover is the intersection point on a ground track where a track ascending in latitude intersects 
a descending track. By differencing altimeter measurements at crossovers, common errors 
cancel leaving the effect of temporal change between the ascending and descending tracks plus 
residual errors in the measurement after corrections. To provide adequate coverage of the West 
Antarctic Ice Sheet, it was determined that the ground tracks should provide coverage to at least. 
86° S. The small uncovered area between 86° S and the South Pole will be completely enclosed 
by the highest density of measurements at the maximum latitude of 86° S. Occasional off-nadir 
pointing to the pole can be performed to provide complete elevation mapping and limited 
elevation change measurements to the pole. Another consideration was the geometry of the 
crossovers, which tend toward larger crossing angles for a 94° retrograde compared to an 86° 
prograde orbit. In the latitude band from 84° to 85° (North or South), the ascending and 
descending tracks will be approximately orthogonal. 


The nominal orbit altitude of 600 km is a compromise between the GLAS instrumental 
preference for lower altitudes and the orbit determination and spacecraft-control preferences for 
higher altitude. For GLAS, the laser power needed to achieve a required signal to noise ratio of 
reflected light is proportional to the fourth power of altitude. On the other hand, a lower altitude 
requires more fuel for more frequent orbit adjustments to compensate for atmospheric drag. 

Also, the accuracy of the precision orbit determinations degrades at lower altitudes due to greater 
perturbations from atmospheric drag and gravity anomalies. The orbit eccentricity will be 
0.0013 with a semi-major axis of about 6970 km. The position of the perigee will remain 
essentially fixed (“frozen orbit”) at the northernmost orbital position. Therefore, the spacecraft 
altitude will be primarily a function of latitude, varying between about 604 km above the 
ellipsoid at 86° N to 622 km above at 86° S. 


The pattern of ground-tracks on the Earth will repeat with a cycle of 183-days (182.758 days and 
2723 revolutions). The separation between ground tracks in the same ascending or descending 
direction will be 14.5 km at the equator, 7 km at 60°, and 2.5 km at 80° latitude. Two sub- 
cycles exist within the 183-day cycle, which produce near-repeat tracks at 8-days and 25-days. 
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9. Surface Elevation Determination 


The measured surface elevation is determined from the altitude of the satellite orbit above the 
Earth minus the range to the surface measured by GLAS, corrected for delays of the laser light in 
the atmosphere. The surface elevation is also corrected for variations of the Earth’s surface 
caused by ocean, solid, load, and polar tides. The single-pulse error budget for ICESat elevation 
measurements summarized in Table 9. 1 applies to ice, land and ocean/lake applications. For this 
error budget, it is assumed that the laser pointing error is 1.5 arcsec and the surface slope is 1°. 
The GLAS altimeter measurement is a scalar distance measurement between a reference point in 
the instrument and the Earth’s surface. When combined with information derived from the 
precision attitude determination process (PAD), the instrumentation provides an altitude vector 
measurement. The precision orbit determination process (POD) yields the position vector of the 
GLAS reference point with respect to the Earth’s center of mass. By adding the altitude vector to 
the position vector, the surface location with respect to the center of mass is calculated. The 
contributing factors to these vectors are dependent on various error sources, which determine the 
overall accuracy of the determined surface location. 

Table 9. 1 Single-Shot Error Budget for ICESat Elevation Measu rements 


GLAS range measurement precision 

10 cm 

Radial orbit determination (POD) 

5 cm 

Pointing determination (PAD) 

7.5 cm 

Atmospheric delay 

2 cm 

Atmospheric forward scattering 

2 cm 

Other (tides, etc.) 

1 cm 

RSS 

13.8 cm 


The individual contributions shown in Table 9.1 are discussed in the following sections. The 
nature of these errors is a combination of random and systematic. The GLAS range measurement 
precision is characterized as random, but most of the other errors are systematic for successive 
pulses. However, the components of the error budget are mostly uncorrelated with each other, so 
the RSS of the error budget is a reasonable representation of the overall error for a single pulse. 
Furthermore, although some of these errors are correlated on successive pulses, the errors for a 
series of measurements over a given location during the life of the mission tend to be random 
rather than systematic. One exception to this is geographically-correlated errors in the absolute 
height of the orbit calculations. However, these errors can be mostly time-independent if the 
orbits are systematically processed with the same gravity model calculation procedures and 
therefore not important for elevation change measurements. Another exception is the 
atmospheric forward scattering (cf. section 3.1 and section 7) if the scattering clouds have 
significant seasonal or interannual variations. 

9.1 Precision Orbit Determination 
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The precision orbit determination (POD) techniques demonstrated with GPS data from 
TOPEX/POSEEDON [e.g., Schutz, 1994; Yunck, et al., 1994] will be applied to ICESat POD. 
These techniques are based on the solution of the satellite equations of motion, generally referred 
to as a “dynamic” approach. The GPS receiver provides measurements of pseudorange and 
carrier phase from the GPS satellites above the antenna ground plane. With measurements on 
two frequencies, ionosphere delays are removed. The few millimeter precision of the carrier 
phase makes it the preferred measurement for POD, even though the measurement requires 
resolution of an ambiguity. Two common approaches to ambiguity resolution are: 1) formation 
of differenced measurements to remove common parameters, such a clock errors, and 2) solve 
for clock parameters using un-differenced measurements. In an alternate “kinematic” approach 
commonly used with aircraft applications the measurement geometry is used to determine 
position without a high-accuracy solution of the equations of motion. However, special 
experiments to resolve ambiguities for aircraft applications cannot be used for a satellite. 

Prelaunch analysis of expected POD errors show that extant gravity models can be expected to 
produce an error in the radial component of ICESat position that exceeds 10 cm [e.g., Rim, et al., 
1999], For example, EGM-96 covariance matrix predicts the error to be about 15 cm. However, 
experience with TOPEX/POSEEDON showed that the GPS measurements could be used to 
“tune” the gravity field and reduce the error contribution from this source. ICESat simulations 
by Rim, et al. [1996, 1999] show that using ICESat GPS data can reduce the gravity errors 
sufficiently to reduce the radial orbit error below the 5-cm value in the error budget. In addition, 
ongoing efforts to improve the gravity model, combined with results from the dedicated gravity 
missions CHAMP and GRACE, should produce further improvements prior to and after the 
ICESat launch. 

9.2 Precision Attitude Determination 

The precision attitude determination (PAD) process calculates the pointing of the laser beam and 
the position of the laser footprint on the Earth’s surface. Data from the stellar-reference system 
(SRS) described in section 6 will be processed to determine the GLAS laser pointing direction. 
The instrument-star tracker (1ST) mounted on the GLAS optical bench will enable determination 
of the benches orientation in space, with respect to a celestial reference frame (CRF) defined by 
the star catalog used with the star camera images. Since the laser-profiling camera (LPC) 
operates at 40 Hz, an image of the far field pattern of the GLAS laser will be obtained for each 
pulse, so the variation of the laser direction with respect to the optical bench will be determined. 
Since the laser pulse is also imaged in the laser-reference camera (LRC), although at a 10 Hz 
rate, further monitoring of the laser direction with respect to the CRF will occur when a star 
appears in the LRC field of view. Simulations of the instrumentation functions and calculations, 
including application of a Kalman filter to the data, have shown that a 1.5 arcsec determination 
of laser pointing direction can be achieved [Bae and Schutz, 2000]. The PAD processing will 
provide the laser pointing direction in the CRF, along with transformation matrices to generate 
the location of the illuminated laser spot in the terrestrial reference frame (TRF). The calculation 
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of the location of the laser footprint on the surface will use the corrected GLAS range 
measurement and the POD and PAD outputs. 


In addition to the instrumentation provided on the GLAS optical bench, the bus star trackers are 
available as a backup. However, the hinges used to mount the GLAS instrument to the bus can 
be expected to distort in response to thermal variations experienced in orbit. These distortions 
will be monitored through analysis of the 1ST and bus tracker data. A model will be developed 
during the mission to predict the distortions and to improve spacecraft pointing control of the 
laser optical bench with respect to the Earth using the bus star trackers. 

9.3 Atmospheric Delay Correction 

The atmospheric delay is calculated by integrating the refractivity, as a function of pressure, 
temperature, and relative humidity, along the ray path [Owens, 1967], At optical wavelengths 
the nadir atmospheric delay, ND, can be written as a sum of the hydrostatic delay, NHD, which 
is a function of surface pressure, Ps, and a residual wet delay, NWD, which is a function of total 
column precipitable water vapor, PW such that ND = NHD + NWD, NHD = 10 6 ki R P s /Ma g m 
and NWD = 10 6 k 2 * R PW/M w , where ki and ki' are constants from the refractivity equation that 
depend on the wavelength of the laser, R is the universal gas constant, Md and M w are the 
molecular weights for dry air and water, and is the mean value of gravity in the atmosphere 
column [Quinn and Herring, 2000]. The g m is expressed accurately in terms of the height, H, in 
meters, and latitude, 0, of the measurement g m = 9.8062 (1 - 0.00265 cos(20)-3.1xlO' 7 
(0.9H+7300)) ms' 2 [Saastamoinen, 1972], 

For a surface pressure of 1000 mbar, the nadir hydrostatic delay is approximately 2.35 m. The 
precipitable water vapor contribution to delay is small, but highly variable. The conversion is 
approximately 0.076 mm of delay per kg/m 2 of precipitable water vapor. In polar regions, 
precipitable water vapor values are generally much less than 10 kg/m 2 , and delay is less than 
1mm. Tropical values are as much as 50 kg/m 2 , leading to delays of 4 mm. For off-nadir angles, 
H, a simple mapping function, D n =ND/ cos (T)) , relates the total delay, D n , to nadir delay ND 
[Marini, 1972] with an estimated accuracy better than 1 mm for pointing angles up to 20°. 
Surface pressure is calculated using the National Center for Environmental Prediction (NCEP) 1° 
gridded fields of atmospheric variables. A hydrostatic equilibrium model of the atmosphere is 
integrated from the upper atmospheric levels to the estimated height of the laser footprint in 
order to calculate surface pressure at the height of the laser footprint. The NCEP total column 
precipitable water vapor is used without modification, because the precipitable water vapor is 
small, particularly in polar regions. The GLAS error budget assumes less than 20 mm single- 
pulse atmospheric delay RMS error. Comparisons of model estimates with polar Automatic 
Weather Station (AWS) data in 1994 show a root-mean-squares (RMS) error of less than 5 mbar, 
after a mean offset is removed to correct for uncertainties in station heights. This error 
corresponds to delay errors of less than 12 mm. The RMS error in the calculated precipitable 
water vapor is less than 2 kg/m 2 , or 0. 15 mm of delay, as estimated using the technique of GPS 
meteorology [Quinn and Herring, 1999]. 
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9.4 Tidal Corrections 


The basic solid Earth and ocean tide corrections required are essentially the same as those 
developed for radar altimetry missions, such as Seasat and more recently TOPEX-Poseidon. 
These corrections and the associated software are fully documented in the International Earth 
Rotation Service (IERS) technical notes [e.g. McCarthy, 1996]. However, for the ICESat 
mission, ocean-load tides and tides on ice shelves are more important than for previous missions. 
Ocean-load tides are the secondary tidal effects on the Earth’s crust caused by the varying 
loading of ocean tides. Ocean-load tides are largest in the continental margins and near-shore 
regions where GLAS will provide continuous altimetry profiles from land to ocean. The largest 
ice sheet elevation changes are likely to be found near the coast. Yi et al. [2000] calculated that 
vertical displacements of several cm due to ocean-load tides are expected along the Antarctic 
margin, using the SPOTL software of Agnew [1996]. However, the details of ocean tide models 
are most uncertain right near Antarctica, because of the paucity of direct observations at high 
latitudes and the poor knowledge of bathymetry around Antarctica limits the accuracy of hydro- 
dynamical ocean tide models [Le Provost et al, 1995]. The Antarctic ice shelves respond to tides 
in extremely complex ways, because of the dynamic effects of water circulation inside the ocean 
cavity underlying the ice. In most instances the cavity geometry is poorly known, and surface 
data records are very sparse. A notable exception is the Amery Ice Shelf where extensive 
measurements has led to the development of an accurate tide model constrained by precise GPS 
observations at the surface [Padman et al., 2000], ICESat data over ice shelves will need to be 
corrected for tidal variations. In turn, ICESat will provide important new data on tides at high 
latitudes. 


10. Calibration and Validation 

The fundamental objective of calibration/validation (CV) is to validate the science data products 
to ensure that appropriate geophysical interpretations can be drawn from the products. The 
identification and removal of non-geophysical artifacts from the data products is an important 
element of the CV process. The GLAS instrument will undergo extensive pre-launch testing for 
calibration, including experiments to characterize the dependency of the calibration parameters 
on environmental factors such as thermal variations. The post-launch phase of CV will verify the 
pre-launch parameters in flight through a variety of experiments, which are described in the 
following sections. In particular, the post-launch CV will assess the existence of biases in the 
measured altitude and pointing parameters with respect to the pre-launch values. 

10.1 Detailed Measurements over Small Target Areas 

Small target areas will be instrumented to verify the GLAS performance and validate the various 
components in the error budget. The target areas will be on flat surfaces, such as White Sands 
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(NM) and Bonneville Salt Flats (UT), but southern hemisphere sites will be used also. These flat 
surfaces provide albedo and surface characteristics that are similar to ice sheets, but also aid the 
verification of the instrument measured surface elevation in the near-nadir pointing direction. 

The computed surface elevation is obtained from the difference between the position of the 
GLAS in space (determined from POD) and the position of the flat surface with respect to the 
same origin and corrected for atmosphere delays. This computed altitude, when compared to the 
measured altitude, provides an evaluation of the GLAS performance and the existence of an 
altitude bias. 

The error budget given in Table 9.1 essentially applies to the verification of the altitude 
measurement using flat surfaces. With near-nadir pointing to a flat surface, the contribution of 
pointing errors to range is small (e.g., 30 arcsec pointing error contributes less than 1 cm). If all 
corrections to the altitude measurement were known precisely, the existence of an altitude bias 
would be evaluated by using several measurements over the flat surface to average out the 
gaussian noise in the measurement. For example, 40 shots (one second, or 7 km along-track) will 
reduce the 10 cm instrument precision to 1.6 cm (10 cm/V40). Provided that variations in the flat 
surface can be characterized at a comparable level over 7 km and that other corrections can be 
made at a comparable level, the flat surface will enable assessment of a range bias at the few 
centimeter level. 

To verify the performance of the pointing determination, small target areas will be instrumented 
to emulate a ground-based CCD camera. This emulation will be accomplished in two ways: 1) an 
overflying aircraft carrying a CCD camera will image the GLAS spot on the surface with respect 
to ground fiducial markers and 2) ground detectors in a grid pattern will serve as a raster-device. 
Ground geodetic surveys using high accuracy processing techniques will provide the positions of 
the fiducial markers and detectors at the centimeter level in the same reference frame used for 
POD and PAD processing. Both approaches are being designed to support determination of the 
location of the illuminated spot to an accuracy compatible with the 1.5 arcsec (4.5 m) pointing 
knowledge requirement. Laboratory testing of the detectors with energy levels comparable to 
those expected for the illuminated GLAS spot have been successfully accomplished and flight 
tests of the airborne technique is in progress. In addition, with appropriate timing elements 
attached to the detectors, the time tag of the GLAS pulses will be verified to better than 1 
microsec. 

The White Sands or Bonneville target areas will be further instrumented to validate the 
atmospheric delay. In addition, a ground-based lidar will enable independent characterization of 
the atmospheric properties during a GLAS overflight. To the extent possible, a ground-based 
satellite laser ranging system (e.g., McDonald Observatory) will be used to verify the POD 
determination of GLAS position during the target area overflight. 

10.2. Aircraft Underflights 
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Aircraft flights both before and after ICESat launch are planned for detailed mapping of 
segments of orbit tracks over the Greenland ice sheet and nearby sea ice, over western US, over 
the Dry Valleys in Antarctica, and possibly over the desert in Saudi Arabia. Surface elevations 
will be measured along tracks ranging from 50 - 300 km long to an accuracy of 10 cm or better 
by NASA’s Airborne Topographic Mapper (ATM) scanning laser [Krabill et al., 1995]. These 
measurements will be compared with surface elevations and surface roughness derived from 
GLAS data. The surveyed segments will be sufficiently wide (approximately 300 m) to include 
most GLAS footprints when ICESat is pointed to the segments, and will include terrain with a 
range of surface slopes and roughness. The effect of ICESat pointing errors on GLAS-derived 
elevations will be determined by the surveyed surface slopes, and consequently should be 
separable from altitude errors. Thus, in addition to an overall assessment of GLAS performance, 
comparison between GLAS and ATM elevations will provide estimates of pointing and altitude 
bias at different locations around the orbit. 

The prime objective of the Greenland surveys is to validate performance over both ice-sheet and 
sea-ice surfaces. Surveyed orbit segments will include surface types ranging from smooth, flat, 
and nearly horizontal to sloping and undulating ice sheet, and crevassed glaciers and rough sea 
ice. Some segments will pass over the narrow (few kms) Greenland outlet glaciers, with 
associated high albedo contrast between the dark rock margins and the glaciers. Because the ice- 
sheet surface elevation is expected to change with time, the Greenland ATM surveys will be 
made as close to simultaneous as possible (within 1 or 2 days) with ICESat passes. These flights 
are planned for May, 2002. 

Other ATM surveys will be over stable terrain (little or no elevation change over time) that can 
provide “benchmark” surfaces for routine GLAS validation throughout the mission. Selected 
surfaces are arid and largely vegetation free. Before launch surveys will be over several arid 
sites in the western US (and possibly in Saudi Arabia). Each survey will include all the 183-day 
Mapping-Phase orbits (about 16 km apart in western US) that lie within approximately 100km X 
100km, over terrain ranging from nearly flat to rough and sloping. In addition, segments of 
nearby 8-day post-launch phase orbits will be surveyed. An ICESat overflight of a surveyed 
orbit segment will occur at least twice per month. In addition, surveys of a range of vegetated 
sites in the northwest US will be used to assess GLAS performance over forest and woodland. 
GLAS performance will be validated at different locations around the globe, primarily to ensure 


surveys of the ice-free areas in the Dry Valleys, near the US station in McMurdo, Antarctica, in 
collaboration with USGS and NSF. In addition, it may be possible to obtain surveys of parts of 
the Amery Ice Shelf and the Australian Desert in a collaborative effort with the Australian 
Antarctic Division and others. 


10.3. Integrated Residual Analysis 



Comprehensive on-orbit calibration of the performance of the GLAS instrument and verification 
of its data products will use data analysis methods based on the techniques of integrated residual 
analysis as described by Rowlands et al. [2000], Residuals are formed between the observed and 
computed parameters in each of the three data types that constitute laser altimetry: range, 
pointing, and orbit. Analysis proceeds by simultaneous processing and weighting of the various 
data types in order to minimize residuals. In this way, all data types both contribute to the 
solution and are refined as systematic errors (biases) are identified and removed from the 
integrated solution. When GLAS observations fall on reference surfaces where the surface 
elevation is known at the 1 to 10-centimeter level, GLAS range data will be used to refine the 
pointing angle data and the ICESat orbit solution, as well as the laser range and timing biases. 
Reference surfaces include the open ocean and land targets with well-characterized and stable 
topography. Single tracks of laser elevation profiles and cross-overs between two independent 
profiles over the reference surfaces are used in the analyses. Rowlands et al. [2000] describes 
specific forms of residual analysis including Range-Residual Analysis (RRA) and Crossover- 
Residual Analysis (CRA). 

The basis of RRA and CRA is a Bayesian least-squares estimation and batch processing 
approach. A rigorous measurement model is employed in which the various parameters of laser 
altimeter pointing and range data, orbit tracking data, spacecraft force models, and geophysical 
models are simultaneously adjusted to minimize the range residuals. All significant terms that 
affect the 3-D location of the laser footprint on the Earth's surface are measured or modeled and 
described in a system of simultaneous differential equations (one for each laser observation). 

The equations are then solved for hundreds of unknown biases and model parameters. Typically 
thousands of data equations are used. The output of the parameter estimation process is a set of 
attitude, range, and orbit bias terms for calibration of the GLAS instrument and the attitude and 
orbit determinations. The methodology utilizes GEODYN, a set of state-of-the-art orbit 
determination and geodetic parameter estimation software [Pavlis et al, 1999]. The RRA 
technique for laser altimetry was initially developed to improve the pointing bias knowledge in 
the Shuttle Laser Altimeter (SLA) experiment [Garvin et al., 1998], The SLA measurements had 
an initial attitude uncertainty from 0. 1° to 1°, because the laser pointing angle was not directly 
measured. The SLA RRA resulted in as much as a two-order of magnitude improvement in 
attitude knowledge [Luthcke et al, 2000a], Continued refinement of the relatively imprecise 
SLA data sets using RRA and CRA is currently yielding range residuals over the ocean surface 
at the 12 cm (1 -ct) level [S. Luthcke, Pers. Com.]. 

For ICESat, roll and pitch maneuvers of the spacecraft, called ocean sweeps, will be periodically 
performed to point the laser beam up to 5° off-nadir to specifically achieve sub-arc sec 
calibration of laser beam pointing angle biases and thus address the largest error term in GLAS 
altimetry. This calibration process employs the rather large sensitivity of laser range 
measurement to laser-beam pointing-angle [Bufton, 1989; Gardner, 1992; Harding et al., 1994] 
and uses range-residual analysis (RRA) with the ocean surface as reference. Open-ocean surface 
elevations are known from radar altimeter measurements and tidal modeling to an accuracy of 3- 
to-12 cm (1 -ct) [Luthcke et al., 2000b], providing a global reference surface throughout the 
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mission. Ocean waves are not a factor due to averaging over thousands of individual laser pulses 
along 100 km to 10,000 km ocean segments. Luthcke et al. [2000b] have shown that the GLAS 
laser attitude can be determined to the sub-arc sec ( 1-a) level from a 30 minute duration (i.e. 1/3 
orbit) ocean sweeps, in which the spacecraft points off nadir in a slow, conical-shaped roll and 
pitch maneuver up to a maximum of 5°. This simulation employed realistic, worst-case error 
models and included the effects of ocean surface waves. For practical considerations, ICESat’s 
ocean sweeps will be diamond-shaped, providing most of the simulated sensitivity. Present 
plans call for one ocean sweep per week minimum, with more frequent sweeps during the initial 
CV period. Biases will thus be assessed as a function of temporally varying parameters, such as 
the orientation of the spacecraft with respect to the sun, which may thermally induce pointing 
variations. The ocean-sweep simulations produce a sub-5cm range bias solution simultaneously 
with the sub-arcsec level attitude bias determination [Luthcke et al., 2000b]. Applying the same 
measurement model as used for ocean sweeps, range-residuals will also be evaluated against 
Digital Elevation Models (DEMs) of stable landform surfaces by airborne laser altimeter 
mapping. 

For ICESat, the CRA is expected to provide cm-level verification of elevation differences at 
cross-overs and improve the estimates of range and attitude biases. As with RRA, CRA uses the 
GEODYN software suite. Constraint equations are formulated in terms of a minimum distance 
between two intersecting altimeter profiles, instead of just the height difference at a crossover 
point. The CRA approach was first developed to exploit the small footprint capability of laser 
altimetry and to take into account the simultaneous improvement in range, attitude, and orbit 
biases. Rowlands et al., [1999] demonstrated an order-of-magnitude improvement in the 
accuracy of Mars Orbiting Laser Altimeter (MOLA) geolocation using CRA. 

Residuals can also be formed between the geolocated laser footprints and the 3-dimensional 
coordinates of land topography at reference sites having a high-resolution DEM. By differencing 
the GLAS elevation profiles with DEM data, the GLAS accuracy can be assessed using methods 
described in Rowlands et al. [2000]. If the DEM has meter-level spatial resolution and 
significant vertical structure within a GLAS footprint, it is further possible to conduct integrated 
residual analysis for pulse waveforms. The time history of the laser backscatter from sloping 
and/or vegetated terrain is compared with a synthesized waveform for validation of footprint 
location and direct comparison with derived parameters (e.g. surface slope and vegetation 
height). This process of waveform matching is described in Rowlands et al. [ 2000] and Blair & 
Hofton [1999], 


11. Data Products and Availability 

The ICESat/GLAS standard data products are summarized in Table 11.1, and more detailed 
information is given at http://glas.wff.nasa.gov/ . The 16 standard data products (GLA01 to 
GLA16) will be created by the ICESat Science Investigator-led Processing System ( I-SIPS) at 
the Goddard Space Flight Center. After the initial calibration/validation period, the data will be 
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processed and products sent within a few weeks following data collection to the National Snow 
and Ice Data Center (NSIDC) ECS DAAC. NSIDC will distribute data requested by the 
community of users. All products are time-ordered, as collected along track, and will be 
distributed in granules (files) as defined in the table. Metadata (time, location, quality, etc) are 
provided with each granule, so the metadata can be searched using NSIDC tools to determine 
which granules to order. Several times during the mission, the data will be reprocessed to 
provide improved products with improved algorithms. 

Table 11.1 ICESat/GLAS Standard Data Products 


Product 

ID 

Granule size 
Mbytes 

Product Name 
GLAS/ICESat-Level..) 

Primary Intended Users 

GLA00 

6 hrs 

L0 Global Instrument 
Packet 

Instrument team, ISIPS 

GLA01 

- l /4 orbit 34 MB 

LI A Global Altimetry 

ISIPS, Science team, 
altimeter scientists 

GLA02 

2 orbits 703 MB 

LI A Global Atmosphere 

ISIPS, Science team 

GLA03 

2 orbits 13 MB 

LI A Global Engineering 

ISIPS, Instrument team 

GLA04 

6 hrs 610 MB 

LI A Global Stellar 
Reference and GPS 

Science and Instrument teams 

GLA05 

orbit 21 MB 

LIB Global Waveform- 
based Range Corrections 

ISIPS, Science team, 
altimeter scientists 

GLA06 

~ Va orbit 9 MB 

LIB Global Elevations 

ISIPS, Science team, 
altimeter scientists 

GLA07 

2 orbits 828 MB 

LIB Global Backscatter 

ISIPS, Lidar scientists, 
Science team 

GLA08 

14 orbits 2 MB 

L2 Global PBL & Elev. 
Aerosol Layer Heights 

ISIPS, Atmospheric scientists 

GLA09 

14 orbits 81 MB 

L2 Global Cloud Heights 
for Multi-layer Clouds 

ISIPS, Atmospheric scientists 

GLA10 

14 orbits 301 MB 

L2 Global Aerosol 
Vertical Structure 

ISIPS, Atmospheric scientists 

GLA11 

14 orbits 12 MB 

L2 Global Thin Cloud/ 
Aerosol Optical Depths 

ISIPS, Atmospheric scientists 

GLA12 

14 orbits 49 MB 

L2 Polar Ice Sheet 
Altimetry 

Ice scientists 

GLA13 

14 orbits 22 MB 

L2 Sea Ice Altimetry 

Ice scientists 

GLAM 

14 orbits 273 MB 

L2 Global Land Surface 
Altimetry 

Land scientists and 
topographic mapping 

GLA15 

14 orbits 240 MB 

L2 Ocean Altimetry 

Ice and Ocean scientists 

GLA16 

14 orbits TBD MB 

Global L3 Elevation and 
Atmosphere (HDF-EOS) 

Multidisciplinary scientists 


The products are defined as level 0, level 1, level 2, and level 3 according to the EOS standards. 
Level 0 is the raw telemetry data, level 1 includes instrument parameters, and level 2 and 3 have 
geophysical ice, ocean, atmosphere, and land parameters. The level 0 products in the EDOS 
standard Production Data Set (PDS) format will be provided directly from the EOS Data and 
Operations System (EDOS) in 6 hour files. The level 1 and level 2 products are in easy-to-read 
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integer binary format. The GLAS team will supply NSIDC with readers in Fortran 90 and EDL 
for all level 1 and level 2 products for distribution to users. For more general users with software 
capable of reading HDF formats, a level-3 global-coverage product is provided in HDF-EOS 
swath format, containing parameters at various data rates (surface elevations, surface height 
distribution, minimum and maximum elevations, cloud and aerosol layers, and optical depths). 
Additional level 3 and 4 products will be defined to include gridded digital elevation models and 
atmospheric (lidar) backscatter images. 

The levels 1 and 2 products are separated into surface-altimeter and atmosphere-related products. 
Level la products are created from level 0 telemetry by applying instrument calibrations and 
converting the data to science units. Level lb are created from level la products. The level 1 and 
2 products are suited for detailed statistical analysis of along-track surface and atmospheric 
properties, temporal changes, and comparisons with other altimetry missions. Level 1 data will 
be mostly of interest to scientists familiar with altimetry or lidar measurements and the 
corrections that are required to derive the along-track geophysical parameters. The level 2 
products contain geophysical and atmospheric parameters with associated instrument and 
geophysical corrections. 

The level 1 altimetry products (GLA01, 05, and 06) are global products at 0.25 orbits per granule 
with 1 second records, each record with one frame of data for 40 laser pulses. The granules are 
started at +/- 50 deg latitude to separate the polar regions. GLA01 contains all telemetered data 
related to altimetry, including the transmitted and received waveforms. GLA05 is an 
intermediate product with information calculated from the waveforms. GLA06, which is 
equivalent to the geophysical data records distributed for radar altimetry missions, contains 
elevations corrected for tides and atmospheric delays, and surface characteristics within the 
footprint. The level 2 elevation products (GLA12-15) are regional products (ice sheet, sea ice, 
land, ocean) at 14 orbits per granule, starting and stopping at the same demarcation as the level 1 
waveform and elevation products. Each regional product contains surface elevation, elevation 
distribution, reflectance, and information related to surface roughness and slope at the 40/s pulse 
rate. 

The level 1 atmosphere products are at approximately 2 orbits per granule, time-ordered in 1 
second records. GLA02 contains the normalized lidar signal, with fundamental corrections and 
normalizations to the raw data, and an estimate of the height of the first bin in the profile and the 
surface elevation. The level 2 atmospheric products are time -ordered 4-second records that 
contain the specific set of parameters related to cloud layers, aerosol layers, cloud and aerosol 
backscatter and extinction profiles, and optical depth. 

12. Summary 

The ICESat mission has a broad mixture of multidisciplinary scientific objectives in glaciology, 
atmospheric science, climatology, and land processes, even though it carries only a single 
instrument. The GLAS instrument has been designed with a combination of synergistic 
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measurement capabilities, without compromising the principal scientific objective of measuring 
small changes in the surface elevations of the polar ice sheets. The small well-located footprint 
of the laser beam, with a well-defined surface reflection, will significantly improve the accuracy 
and coverage of satellite altimeter measurements compared to prior satellite radar altimetry. 

Also, the dedication of the mission to a primary scientific objective required a high-latitude orbit 
inclination that is optimized for ice sheet research. 

The inclusion of the lidar channel in GLAS required approximately 50% higher power in the 
primary 1064 nm laser channel and additional detectors and electronics, while utilizing the same 
telescope and many other common elements of the instrument, spacecraft, and the overall 
mission. In return, the GLAS atmospheric measurements not only broaden the science of the 
ICESat mission to studies of the climatic impacts of clouds and aerosols, but will also improve 
out understanding of the effects of the atmosphere on the accuracy of the ice sheet altimetry 
measurements. In addition, ICESat’s atmospheric measurements should greatly improve our 
knowledge of the distribution, heights, and types of polar clouds. Previous measurements of 
cloud cover in polar regions have been very limited by the capabilities of spacebome passive 
sensors in the polar darkness. The inadequacy of such knowledge, for example, seriously limited 
our ability to make accurate estimates of the percent of polar cloud cover than will obscure the 
GLAS primary altimeter range measurements over the ice sheets. Also, of direct importance to 
ice sheet science, will be improvements in our understanding of polar clouds and precipitation, 
which is a major determinant of the surface mass balance of the ice sheets. While GLAS could 
have been designed to operate only over the polar ice sheets, global operations enables 
measurements of land and ocean topography for significant improvements in the accuracy of 
land topography and a variety of applications that will develop as various investigators explore 
the utility of ICESat data. 

Initially, the major satellites in NASA’s Earth Observing System (EOS) were each planned as a 
series of three 5-year missions to systematically acquire global-change data sets for a period of 
15 years. The first mission, Terra, was launched in December 1999 and the second mission. 
Aqua, should be launched a few months before ICESat. The schedules and description of these 
and associated missions are given at http://www.earth.nasa.gov/ along with a research strategy 
for follow-on missions. One category of observations is systematic measurements that are 
needed for long-term continuity of data on critical environmental variables inherent in global 
change. The NASA’s research strategy notes that systematic is not necessarily synonymous with 
continuous measurement, but continuous measurements are required when short-term natural 
variability or calibration uncertainties between successive discontinuous records would mask 
significant long-term trends. For example, overlapping measurement records from successive 
sensors are required when no ground-based observation can provide an independent calibration 
standard. Although the ERS- 1 and ERS-2 missions and their radar altimeters are nearly 
identical, their measurements over ice sheets were found to have a spatially-variant inter-satellite 
bias of 5 to 40 cm that did not occur over oceans and could not have been discovered if the 
mission operations had not overlapped. For ICESat, the calibration plan is intended to reduce the 
vulnerability to gaps in successive laser altimetry mission. More important, however, are the 
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seasonal and interannual variability in ice sheet surface elevations and mass balance. For 
examples, the data in Figures 2.2 for 7 years of elevation change and Figure 2.4 for 20 years of 
temperature and surface melting illustrate a well-know inter-annual variability in ice-sheet 
surface mass balance and surface climate that might preclude deductions on decadal-scale trends 
from records of short duration or with gaps of several years. 
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Figure 1 . 1 ICESat will carry the Geoscience Laser 
Altimeter System (GLAS) for measurement of 
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Figure 2.3. Elevation change map of a portion of the West Antarctic 
ice sheet showing up to -30 cm/yr thinning in the Pine Island/Th waites 
drainage basin in the Northeast and up to 25 cm/yr thickening along the 
ridge between the Pine Island/Thwaites Glacier and the Ross Ice-Shelf 
drainage basin in the western part. dH/dt values mapped on 50 km 
grid are from multi-variate linear and sinusoidal fits to H(t) time series 
of elevation differences at crossovers within 100 km radius and ± 250 
m elevation of the grid center. 
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Figure 2.4. Mean extent of surface melting on the Greenland ice sheet 
and average temperature anomalies at six coastal stations for 
1979-1999. The increasing trend in melt extent is 0.97%/yr while the 
temperatures warmed 0.24 C over the 21 -year period. (Melt- 
temperature correlation coefficient is 0.82.) Melt measurements are 
from long-term satellite passive- microwave sensing. 


6 






0 1 000 2000 3000 

Distance (km^ 


Figure 3. 1 A nighttime (upper) and daytime (lower) simulation of the GLAS 532 nm 
signal using aircraft lidar data acquired during the TOGA-COARE experiment. 

Color coding represents the expected received signal in photons per height bin (76.8 
meters/bin). 
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Figure 6. 1 - Views of GLAS showing at top the 3 laser 
boxes (yellow) on the optical bench, 1 - m diameter 
telescope with shroud, heat pipe(red), and side radiators, 
and at bottom the GLAS star tracker (pink), electronics 
boxes, and the small telescope (grey) of the stellar 
reference system. 
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Ranging Error vs. Received Signal Photons/Pulse 


Laser pulsewidth: 5ns FWHM 

Waveform digitizer: IGs/s, 8 bits, signal 1/4 full scale 

March 18, 1999 



Figure 6.2. Theoretical estimate of the GLAS ranging errors as a function of the detected photons/pulse for < 
ratios and various pulse spreadings from surface slopes. The expected signals for clear-sky atmospheric tran 
photons/pulse from terrain and 20,000 photons/pulse from snow/ice surfaces. 
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Figure 7.1. Calculation of range and elevation distribution from typical transmitted and 
received GLAS pulses. 
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Laser Pulse Broadening toy Surface Type 
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Figure 7.2. Characteristics of returned laser pulse as a function of surface slope and 
roughness, both of which broaden the pulse. 



